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With a sample of approximately 89 million BB pairs collected with the BAS4R detector, we 
perform a search for B meson decays into pairs of charmless vector mesons (0, p, and K*). We 
measure the branching fractions, determine the degree of longitudinal polarization, and search for CP 
violation asymmetries in the processes S"*" — » (j)K*~^ , B° — > (pK*'^, B^ — > p'^K*^ , and B'^ — > p^p"*". 
We also set an upper limit on the branching fraction for the decay S° — > p^p". 
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Charmless B meson decays provide an opportunity 
to measure the weak-interaction phases arising from 
the elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
quark-mixing matrix and to search for phenomena 
outside the standard model, including charged Higgs 
bosons and supersymmetric particles ^ . 

The decays to two vector particles are of special inter- 
est because their angular distributions reflect both strong 
and weak interaction dynamics 0- The asymmetries 
constructed from the number of B decays with each fla- 
vor and with each sign of a triple product are sensitive to 
CP violation or to final state interactions (FSI) 0. The 
triple product is defined as (qi — q2) ■ Pi x p2, where qi 
and q2 are the momenta of the two vector particles in 
the B frame and pi and P2 represent their polarization 
vectors. 

The first evidence for the decays of B mesons to pairs 
of charmless vector mesons was provided by the CLEO 5] 
and BABAR experiments with the observation of i? 
(j)K* decays. The CLEO experiment also set upper limits 
on the B decay rates for several other vector- vector final 
states 0. The BELLE experiment recently announced 
observation of i3+ p^p^ 

In this analysis, we use the data collected with the 
BABAR detector |9| at the PEP-II asymmetric-energy 
e+e" collider These data represent an integrated 

luminosity of 81.9 fb^^, corresponding to 88.9 million 
BB pairs, at the T{AS) resonance (on-resonance) and 
9.6 fb~^ approximately 40 MeV below this energy (off- 
resonance). The T{AS) resonance occurs at the e+e" 
center-of-mass (cm.) energy, y^, of 10.58 GeV. 

Charged-particle momenta are measured in a tracking 
system that is a combination of a silicon vertex tracker 
(SVT) consisting of five double-sided detectors and a 40- 
layer central drift chamber (DCH), both operating in a 
1.5 T solenoidal magnetic field. Charged-particle iden- 
tification is provided by the energy loss (dE/dx) in the 
tracking devices (SVT and DCH) and by an internally re- 
flecting ring-imaging Cherenkov detector (DIRC) cover- 
ing the central region. Photons are detected by a CsI(Tl) 
electromagnetic calorimeter. 

We search for charmless vector- vector B meson decays 
involving (p, p, and if* (892) resonances. The event se- 
lection and analysis technique have been discussed ear- 
lier 0. We fully reconstruct the charged and neutral 
decay products including the intermediate states — > 
K+K', K*° K+TT- and K°Tr°, K*+ K+w" and 
K'^tt^, TT+TT^, — > TT^TT*^, with TT*' — !■ 77 and 

—)■ Kg —^ vr+TT^, where inclusion of the charge conju- 
gate states is implied. Candidate charged tracks are re- 
quired to originate from the interaction point. Looser cri- 
teria are applied to tracks forming Kg candidates, which 
are required to satisfy |m^+^- — mj^o\ < 12 MeV with 
the cosine of the angle between their reconstructed flight 
and momentum directions greater than 0.995, and the 
measured proper decay time greater than 5 times its un- 



certainty. Charged-particle identification provides sepa- 
ration of kaon tracks from pion and proton tracks. 

We reconstruct tt*^ mesons from pairs of photons, each 
with a minimum energy of 30 MeV. The invariant mass 
of the 7r° candidates is required to be within 15 MeV of 
the nominal mass. The helicity angle of a (f>, K* , or p is 
defined as the angle between the momentum (pi or P2) of 
one of its two daughters (if"*", K, or tt+j respectively) in 
the resonance rest frame and the momentum (qi or q2) of 
the resonance in the B frame. To suppress combinatorial 
background with low-energy tt" candidates we restrict 
the K* — > Ktt^ and p+ — > tt+tt'^ helicity angle range to 
cose* < +0.5. 

We identify B meson candidates kinematically using 
two nearly independent variables 9]: the bcam-cnergy- 
substituted mass toes = [(s/2 + p^ • PB)^/Ef — p|]^/^ 
and the energy difference AE = [EiEs — PrPs — 
s/2)/y/s, where {Ei,pi) is the initial state four- 
momentum obtained from the beam momenta, and 
(Eb ,Pb) is the four-momentum of the reconstructed B 
candidate. Our initial selection requires toes > 5.2 GeV 
and |A£;| < 0.2 GeV. 

To reject the dominant quark-antiquark continuum 
background, we require | cos 61^1 < 0.8, where 6t is the 
angle between the i?-candidate thrust axis and that of 
the rest of the tracks and neutral clusters in the event, 
calculated in the cm. frame. We also construct a Fisher 
discriminant that combines eleven event-shape variables 
defined in the cm. frame [^lll|. 

Monte Carlo (MC) simulation ^ demonstrates that 
contamination from other B decays is negligible for the 
modes with a narrow (j) resonance and is relatively small 
for other charmless B decay modes. We achieve fur- 
ther suppression of i3-decay background by removing sig- 
nal candidates that have decay products consistent with 
D — !■ Kir, Kirn decays. The remaining small background 
coming from B decays (about 6% of the total back- 
ground) is taken into account in the fit described below. 
In this analysis we do not explicitly provide a fit com- 
ponent for other partial waves with the same final-state 
particles selected within vector resonance mass windows. 

We use an unbinned, extended maximum-likelihood 
(ML) fit to extract signal yields, asymmetries, and 
angular polarizations simultaneously. We define the 
likelihood Ci for each event candidate i as the sum 
of njkPj{xi]a) over three event categories j, where 
Vj{xi\a) are the probability density functions (p.d.f.) 
for measured variables Xi, rijk are the yields to be ex- 
tracted from the fit, and k is the measured tag (1 or 2, as 
defined for asymmetry measurements later). There are 
three categories: signal {j = 1), continuum qq {j = 2), 
and BB combinatorial background {j = 3). The fixed 
numbers a parameterize the expected distributions of 
measured variables in each category. They are extracted 
from MC simulation, on-resonance Ai?-TOEs sidebands, 
and off-resonance data. 
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TABLE I: Summary of results for the measured _B-decay modes; e denotes the reconstruction efficiency and etot the total 
efficiency including daughter branching fractions, risig is the fitted number of signal events, B is the branching fraction, 
is the longitudinal polarization, and Acp is the signal charge asymmetry. The decay channels of K* are shown when more 
than one final state is measured for the same B decay mode. All results include systematic errors, which are quoted following 
the statistical errors. The errors are combined for the reconstruction efficiency. The upper limit on the — > p^p" branching 
fraction is given at 90% confidence level including systematic uncertainties and conservatively assuming the efficiency for /l = 1. 



Mode 


e (%) 


etot (%) 




B (xlO-6) 


II 


Acp 


->kK+ 


23.9 ±2.1 
14.3 ±1.4 


5.0 

2.7 
2.3 


33.3t6;4 ± 1.2 
22.3tg;5 ± 3.2 


12.7 l2;o =t 1-1 

13.9 ± 1.2 
10.7 ^3;^ ± 1.8 


0.46 ±0.12 ±0.03 
0.50 ±0.03 
0.40 lo;ig ±0.06 


+0.16 ±0.17 ±0.03 
-0.02 ±0.20 ±0.03 
+0.63 lo;3i ±0.05 




29.7 ±2.6 
10.5 ± 1.0 


10.3 
9.7 
0.6 


101 t]i ± 3 

2.ot5;3 ± 0-6 


11.2 ± 1.3 ± 0.8 
11.7 ±1.4 ±0.8 
3.8 


0.65 ±0.07 ±0.02 
0.64 ±0.07 ±0.02 
1.00 to.se =•= 0.25 


+0.04 ±0.12 ±0.02 
+0.04 ±0.12 ±0.02 




12.3 ±2.0 
6.0 ± 1.4 


4.8 
2.8 
2.0 


35.7tn;o =1=3.6 


lO.eljIg ± 2.4 
14.3l4;4 =•= 2-9 
4.8l^:?±3.2 


0.96 toll ± 0.04 
0.90 to.le =•= 0-04 
1.00 to.2o =1=0.03 


+0.20 to 29 ± 0.04 
+0.17 ±0.04 
+0.28 to g2 ±0.19 




4.7 ±0.9 


4.6 


93 =•= 10 


22.5tg;4 ±5.8 


0.97 l°;o7 ±0.04 


-0.19 ±0.23 ±0.03 




17.6 ± 1.5 


17.6 


9.7t}li ± 2.0 


< 2.1 







The fit input variables Ai?, mEs, Fisher dis- 

criminant, invariant masses of the candidate K* and (p 
(or p) resonances, and the K* and (j) (or p) helicity angles 
9i and 62 ■ The correlations among the fit input variables 
in the data and signal MC are found to be small (typ- 
ically less than 5%), except for angular correlations in 
the signal as discussed below. The p.d.f. Vj{xi]a) for a 
given candidate i is the product of the p.d.f. 's for each 
of the variables and a joint p.d.f. for the helicity angles, 
which accounts for the angular correlations in the signal 
and for detector acceptance effects. We integrate over the 
angle between the decay planes of the two vector-particle 
decays, leaving a p.d.f. that depends only on the two he- 
licity angles and the unknown longitudinal polarization 
fraction = Fi/F. The differential decay width 3j 
d^F/d cos Videos 02 is 

^ I ^(1 - /l) sin' 0^ sin^ 62 + /l cos^ 6^ cos^ • (1) 

To describe the signal distributions, we use Gaussian 
functions for the parameterization of the p.d.f. 's for /S.E 
and mESj and a relativistic P-wave Breit-Wigner distri- 
bution, convoluted with a Gaussian resolution function, 
for the resonance masses. For the background, we use 
low-degree polynomials or, in the case of mssj £^n empir- 
ical phase-space function [l^j- The background parame- 
terizations for resonance masses also include a resonant 
component to account for resonance production in the 



continuum. The background helicity-angle distribution 
is again separated into contributions from combinatorial 
background and from real vector mesons, both described 
by polynomials. The p.d.f. for the Fisher discriminant 
is represented by a Gaussian distribution with different 
widths above and below the mean. 

We denote Qtp as the sign of the triple product and 
Qch as the B- flavor sign (Qch = +1 for B and Qch = — 1 
for B). The charged B is intrinsically flavor-tagged. The 
flavor of a neutral B is determined from the charge of the 
kaon in the final states with the K*^ K^tt~, but is 
undetermined for the decay mode K*'^ X^tt" and for 
the decay B° p°p°. 

We rewrite the event yields rijk (fc=l,2) in each cate- 
gory j in terms of the asymmetry Aj and the total event 
yield Uj: Uji = nj x {l+Aj)/2 and nj2 = nj x {l — Aj)/2. 
We define three signal asymmetries using the tags k: Acp 
{k = 1 for Qch > 0, fc = 2 for Qch < 0), Ap {k = 1 for 
Qch X Qtp > 0, fc = 2 for Qch x Qtp < 0), and Ap 
{k — 1 for Qtp > 0, fc = 2 for Qtp < 0). A non-zero 
value for Acp would provide evidence for direct- CP vi- 
olation, non-zero Atp indicates CP violation even in the 
absence of FSI, and y^gp is sensitive to strong-interaction 
phases Q. 

We allow for multiple candidates in a given event by as- 
signing to each a weight of 1/iVi, where Ni is the number 
of candidates in the same event. The extended likelihood 
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for a sample of A''cand candidates is 



exp 




(2) 



The event yields rij , asymmetries Aj , and polarization 
/l are obtained by minimizing the quantity = — 2 ln£. 
The dependence of on a fit parameter rij, Aj, or 
is obtained with the other fit parameters floating, their 
values are constrained to the physical range < /l < 1 



and < 



We quote statistical errors corresponding 



to unit change in x^- When more than one K* decay 
channel is measured for the same B decay, the channels 
are combined by adding their distributions for rij , Aj , 
or /i. The statistical significance of a signal is defined as 
the square root of the change in when constraining the 
number of signal events to zero in the likelihood fit. If 
no significant event yield is observed, we quote an upper 
limit for the branching fraction obtained by integrating 
the normalized likelihood distribution. Performance of 
the ML fit is tested with generated MC and control sam- 
ples. 

The results of our maximum likelihood fits are sum- 
marized in Table ^ For the branching fractions, we as- 
sume equal production rates of B'^S'^ and B^B^ . We find 
significant signals in p^K*~^ (4.8ct), (6.1a), and in 

both (j)K* (above lOcr each) decay modes. We measure 
the charge asymmetries and longitudinal polarizations in 
the above modes. The projections of the fit results are 
shown in Fig. ^ find [3 The asymmetries involving triple 
products are obtained from separate fits. The results are 
shown in Table ITTl 

Systematic uncertainties in the ML fit originate from 
assumptions about the p.d.f.'s. We vary the p.d.f. pa- 
rameters within their respective uncertainties, and derive 
the associated systematic errors. The signals remain sta- 
tistically significant under these variations. Additional 
systematic errors in the number of signal events origi- 
nate from uncertainty in the background component for 
pK* that peaks in m^s, where we take the uncertainties 
to be the estimated values. 

The systematic errors in the efficiencies are for track 



TABLE IL Summary of asymmetry results with triple prod- 
ucts discussed in the text. 



Mode 



Atp 



Asp 



-0.02 ±0.18 ±0.03 
+0.06 ±0.12 ±0.02 
+0.03 ±0.29 ±0.03 
+0.09 ±0.24 ±0.04 



-0.04 ±0.18 ±0.03 
+0.07 ±0.12 ±0.02 
+0.28 t°;33 ±0.04 
-0.23 ±0.24 ±0.04 
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FIG. 1: Projections of the multidimensional fit onto the 
variable tues for (a) B+ -> (b) B° -» (j)K*°, (c) 

p"K'"^ , and (d) —^ p"p^ candidates after a require- 
ment on the signal-to-background probability ratio "Psig/Pbiig 
with the p.d.f. for jties excluded. The points with error bars 
show the data, the solid (dashed) line shows the signal-plus- 
background (background only) p.d.f. projection. 
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FIG. 2: Invariant mass projections (a) 0, (b) K* for B 
(l>K'- (c) p°, (d) K'-+ for B+ -> p°K*+- (e) p°, (f) p+ for 
B^ p^p^ candidates after a requirement on the signal- 
to-background probability ratio Psig/Pbitg with the p.d.f. for 
mass excluded. For point and line definitions see Fig.0 



finding (0.8% per track), particle identification (2% per 
track), and Kg and tt" reconstruction (5% each). Other 
minor systematic effects are from event-selection crite- 
ria, daughter branching fractions statistics, and 
number of B mesons. We account for the fake com- 
binations in signal events passing the selection criteria 
with a systematic uncertainty of 3-12%, depending on 
the mode. The reconstruction efficiency depends on the 
decay polarization. We calculate the efficiencies using 
the polarization measured in each decay mode |l5l| (com- 
bined for the two (j)K* modes) and assign a systematic 
error corresponding to the total polarization measure- 
ment error. The B'^ — > p'^p'^ branching fraction limit 
incorporates uncertainties in the p.d.f.'s and in the re- 
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construction efficiency, while we conservatively assume 
/l = 1 for the efficiency (which is 29% for /l = and 
18% for A = 1). 

In the polarization and asymmetry measurements, we 
again include systematic errors from p.d.f. variations 
that account for uncertainties in the detector acceptance 
and background parameterizations. The biases from the 
finite resolution in helicity angle measurement and di- 
lution due to the presence of the fake combinations are 
studied with MC simulation and are accounted for with 
a systematic error of 0.02 for polarization. We find the 
uncertainty on the charge asymmetry due to the track 
reconstruction efficiency to be less than 0.02 0. The 
asymmetry measurements are corrected by the small di- 
lution factors. 

In summary, we have observed the decays — )■ 
(j)K*+, B° (j)K*^, B+ p°K*+, and B+ p^p+, 
measured their branching fractions and longitudinal po- 
larizations, and looked for asymmetries sensitive to CP 
violation and FSI. These results supersede the earlier 
BABAR measurements of the B — > (pK* . Our asymme- 
try results rule out a significant part of the physical re- 
gion, providing constraints on models with hypothetical 
particles, but are not yet sufficiently precise to allow de- 
tailed comparison with standard model predictions. Our 
measurement of longitudinal polarization is of interest for 
the study of decay dynamics. 
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